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ABSTRACT 

In a recent paper the authors presented an extended study aiming at simulating the 
classical meteorological parameters and the optical turbulence at Dome C during the 
winter with the atmospherical mesoscale model Meso-NH. The goal of that paper was 
to validate the model above Dome C with the support of measurements and to use it 
afterwards above the Internal Antarctic Plateau to discriminate between the qualities 
of different potential astronomical sites on the plateau. A statistical analysis has been 
presented and the conclusions of that paper have been very promising. Wind speed 
and temperature fields (important for the computations of the optical turbulence pa- 
rameters) revealed to be very well reconstructed by the Meso-NH model with better 
performances than what has been achieved with the European Centre for Medium- 
Range Weather Forecast (ECMWF) global model, especially near the surface. All 
results revealed to be resolution-dependent and it has been proved that a grid-nesting 
configuration (3 domains) with a high horizontal resolution (AX = 1km) for the in- 
nermost domain is necessary to reconstruct all the optical turbulence features with a 
good correlation to measurements. High resolution simulations provided an averaged 
surface layer thickness just ^14 m higher than what is estimated by measurements, 
the seeing in the free atmosphere showed a dispersion from the observed one of just a 
few hundredths of an arcsecond (Ae ~ 0.05"). The unique limitation of the previous 
study was that the optical turbulence in the surface layer appeared overestimated by 
the model in both low and high resolution modes. In this study we present the results 
obtained with an improved numerical configuration. The same 15 nights have been 
simulated, and we show that the model results now match almost perfectly the obser- 
vations in all their features: the surface thickness, the seeing in the free atmosphere as 
well as in the surface layer. This result permits us to investigate now other antarctic 
sites using a robust numerical model well adapted to the extreme polar conditions 
(Meso-NH). 
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1 INTRODUCTION 

The ability of the model to forecast the evolution of the 
atmosphere and the optical turbulence above the Antarctic 
Plateau has been extensively disc ussed in a previous paper 
by our team l|Lascaux et al.l l2009). In that paper all the 15 
winter nights, for which measurements of the C% are avail- 
able (Trinquet et al. 2008) above Dome C, h ave been simu- 
lated with Meso-NH. The main conclusion of lLascaux et all 
(2009) was that Meso-NH was able to reconstruct the opti- 
cal turbulence in a region with extreme meteorological con- 
ditions such as Dome C. However a mono-domain configura- 
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tion with low resolution (100 km) is not suitable to provide 
optical turbulence features well correlated to measurements 
while a grid-nesting configuration done with three domains 
(horizontal resolutions of 25 km, 10 km and 1 km) does it. 
More precisely, the mean simulated surface layer thickness 
h s i,mnh-high='i8-9 m ± 7.6 m was only 14 m higher than 
the observed one (h 3 i t0 b a = 35.3 m ± 5.1 m). The median 
simulated free atmosphere seeing (e mn h,FA=0-35 ± 0.24 arc- 
sec) was very well correlated to the observed one (e b s ,FA = 
0.3 ± 0.2 arcsec). However the model tended to overes- 
timate the intensity of the optical turbulence in the sur- 
face layer, thus generating a too strong median total seeing 
{£mnh,TOT= 2.29 ± 0.38 arcsec) if compared to the observed 
one (e b s .TOT~i-G ± 0.2 arcsec). This paper is intended to 
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Figure 1. Orography of Antarctica as seen by the Meson-Nh model (polar stereographic projection, grid-nesting configuration), (a), 
(b) and (c) show the original orography (GTOPO30) with horizontal resolution of 25 km, 5 km and 1 km, respectively, (d), (e) and (f) 
show the orography obtained with the new digital elevation model (RAMPDEMv2) with horizontal resolution of 25 km, 5 km and f km, 
respectively. The dot labeled 'C indicates the Concordia Station. The dot labeled 'A' indicates Dome A. The altitude is expressed in 
meter (m). 



present the updated results for the same set of winter nights 
with an improved numerical configuration for the Meso-NH 
model. 



2 IMPROVEMENTS IN THE NUMERICAL 
CONFIGURATION 

We refer the readers to lLascaux et al] |2009) for a complete 
overview of the numerical configurations of the mesoscale 
simulations. Here we briefly describe the differences in the 
numerical set-up with respect to the lLascaux" et all <|2009h 
paper. 

(1) A different digital elevation model (DEM), the so-called 

RAMPDEMv2 (Radarsat Antarcti c Mappi ng Project 

Digital Elevation Model, version 2 (|Liu et all [2001)) has 
been used instead of the GT OPO30 DEM fr om th e U.S. 
Geological Survey used in lLascaux et al.l |2009T ). The 
improvement on the description of the orography is shown 
on Fig. [1] The best description of the Dome C area is 
obtained with RAMPDEMv2, and is especially visible at 
horizontal resolutions below 5 km. At AX=5 km (Fig. [Tp,f, 
with GTOPO30 and RAMPDEMv2 orographic models, 
respectively), one can see that the Concordia Station 
(labelled C on the figure) is located on the local summit 
with RAMPDEMv2, and is not with GTOPO30. This is 
even more evident inside the innermost domain (AX=1 km, 
Fig.Eb,g). 

(2) The surface scheme of Meso-NH, the so called ISBA 
(Interaction Soil Biosphere Atmosphere) scheme, has been 
optimized for antarctic conditions: the thermal coefficient 
of the soil has been optimized for polar conditions and 



a climatolo gical underground tem perature T c has been 
introduced (|Le Moigne et all 120091 ) . We briefly summarize 
the main concepts of that study. Our ability in well 
reconstructing the surface temperature T s is related to 
the ability in well reconstructing the sensible heat flux 
H that is responsible for the buoyancy-driven turbulence 
in the surface layer. The Meso-NH surface scheme is 
based on the force-restore method which consists in two 
equations that control the temporal evolution of the surface 
temperature T s and the deep temperature T2 at a few tens 
of centimeters underground. The equation of the temporal 
evolution of the deep temperature has been modified adding 
a term depending on two free parameters: a climatological 
temperature T c and a relaxation term called 7. The two 
equations previously mentioned were forced by a one-year 
set of measurements of solar direct and long-wavelength 
radiations, the temperature T, wind speed V, pressure p 
and specific humidity q of the air above the ground. The 
simulated T 3 and T2 were compared to the observed one 
(i.e. the temperature measured at -5 cm and -30 cm from 
the ground) up to minimize the dispersion. This permitted 
to determine the two free parameters T c and 7. Such a 
study permitted us to optimize the surface scheme ISBA 
for applications of the Meso-NH model to polar conditions. 

Coming back to our study, both mono-domain (low 
resolution) and grid-ne sting (high reso l ution) configurations 
have been employed in lLascaux et al.1 {2009). In this paper 
we could have restricted the study to the grid-nesting 
configuration (with high horizontal resolution) as it was the 
one giving the best results. However we decided to present 
also the results of the low horizontal resolution to confirm 
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the previous conclusions of lLascaux et all (|2009f ) regarding 
the horizontal resolution necessary to be used for obtaining 
reliable forecasts. 



3 OPTICAL TURBULENCE FORECAST 

We focus our attention on the prediction of the 3D maps 
of C% obtained with the Astro-Meso-NH package imple- 
mented by our team in the Meso-NH mesoscale model (Mas- 
ciadri et al. (1999a), (1999b)) and validated in successive 
phases above different astronomical sites (among the most 
important Masciadri & Jabouille (2001), Masciadri et al. 
(2004), Masciadri & Egner (2006)). The model is employed 
here with the new configuration mentioned in the previous 
section. We especially look at three different features that 
m ainly characterize the optical turbulence, already discussed 
in lLascaux et al" | (I2009T ): 

• the surface layer thickness; 

• the median free atmosphere seeing; 

• the median total seeing (from the ground up to 10 km 
above ground level). 

The st udy is performed for t he same set of 15 winter nights 
used in lLascaux et al.l (|2009T ) and employing the same anal- 
ysis. We first present the results of prediction of the surface 
layer thickness and then we present the forecasting of the 
free atmosphere seeing efa and the total seeing etot- 



3.1 Surface layer thickness (h s i) 

The same method used in lLascaux et al.l (|2009h is employed 
to determine the surface layer thickness of each nights. We 
highlight that this is not the only way to define the sur- 
fac e layer but we are fo rced to use this method employed 
bv lTrinauet et all |2008) to be able to compare simulations 
with measurements. The thickness h s i is defined as the verti- 
cal slab containing 90% of the optical turbulence developed 
in the first kilometer: 



i 8 h : ciwdh 



XI km „. 
3r, C 



%{h)dh 



< 0.90 



(1) 



where C N is the refractive index structure parameter. 

We report in Table [1] the values of the observed sur- 
face layer thicknesses and the corresponding mean and stan- 
dard deviation (a) and statistical error (a/VW) for the 15 
nights. The simulated h s i for both mono-domain configu- 
ration (i.e. low horizontal resolution) and grid-nesting con- 
figuration (i.e. high horizontal resolution) are summarized 
on Table [5] The values related to each night as well as the 
mean h a i, mnh - low and h s ^ mnh - high are computed. A careful 
examination of the night of September 5, 2005 (for which 
the computed h a i tmrl h-high was 338.4 m) highlighted the 
fact that the value of the surface layer retrieved with the 
criterion reported in Eq[T]was unexploitable for this night. 
Indeed a look at the temporal evolution of the C% profile 
during that night (not shown here) clearly indicates a very 
weak turbulent night, with maxima of C% concentrated well 
below 30 m. Using the criterion expressed by Eq[T]it follows 
that h s i has to be necessarily large to achieve the 90% of 
the turbulence developed in the first kilometer in this night. 



Table 1. Surface layer thickness h s [ for 15 winter nights (Trin- 
quet et al. 2008). The mean value is reported with the associated 
standard deviation (a) and statistical error (cr/v^V) where N is 
the number of independent estimates that is independent nights. 
Units in meter (m). 



Date 


Observed surface 


Date 


Observed surface 




layer thickness 




layer thickness 


04/07/05 


30 


12/08/05 


22 


07/07/05 


21 


29/08/05 


47 


11/07/05 


98 


02/09/05 


41 


18/07/05 


26 


05/09/05 


20 


21/07/05 


47 


07/09/05 


39 


25/07/05 


22 


16/09/05 


24 


01/08/05 


40 


21/09/05 


22 


08/08/05 


30 








Mean 


35.3 






a 


19.9 








5.1 
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Figure 2. Correlation plot between measured and simulated sur- 
face layer thicknesses (black: mono-domain configuration; red: 
grid- nested configuration). For the simulated values only the 
mean values between 12 UTC and 16 UTC are considered. For 
each configuration of the simulation (high and low horizontal res- 
olution) the error bars are reported for one point only (and are 
equal to cr). Units are in meter (m). 



That's why we decided to remove this night from the mean 
computation in the high resolution case in order to have 
a more representative h a i tmn h-hi g h- Anyway, this removal 
does not affect the qualitative conclusions inferred from the 
comparison between the observed and the computed mean 
values. 

The high resolution configuration gives a mean surface 
layer thickness h a i tmn h-hi g h = 44.2 ± 6.6 m (with a = 24.6 
m) and the low resolution configuration a h s i : , nn h-iow = 
68.4 ± 6.8 m (with a = 26.5 m). We confirm, ther efore, the 
same conclusion obtained in lLascaux et all (|2009l ): the sur- 
face layer thickness computed in grid-nested mode is closer 
to the observed one (/i s ;, D i>s — 35.3 ± 5.1m). Figure [2] shows 
the correlation between observed and simulated surface layer 
thicknesses with the associated intrinsic dispersion (error 
bar) equal to the a values. This figure confirms the tendency 
of the high resolution simulations to give h s i values closer to 
the observed ones from a statistical point of view. One can 
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Table 2. Mean surface layer thickness h s i for the same 15 win- 
ter nights that are reported in Table [l] deduced from Meso-NH 
computations using the criterion in Eq. [T] in the temporal range 
12:00-16:00 UTC. The mean value is reported with the associated 
statistical error a/\/N. Units in meter (m). 



Date 


Surface layer thickness 


Surface layer thickness 




l\/l OCA l\T 1—1 fTT"1/"l vi not - 1 VI tr 
1V1L.6U IN 11 H,I 11LB L1J.1H, 


A/1 flcn i\T T-i in nun nwi oin 

IVlCBU IN 11 lIlUllU v_i.UlIlell.il 




OU.'l 


JU.O 


07/07/05 


35.4 


37.2 


11/07/05 


80.0 


108.3 


18/07/05 


49.7 


82.2 


21/07/05 


66.7 


94.9 


25/07/05 


27.4 


36.7 


01/08/05 


22.6 


39.6 


08/08/05 


34.2 


71.9 


12/08/05 


16.7 


43.9 


29/08/05 


91.4 


106.1 


02/09/05 


70.9 


89.5 


05/09/05 


338.4 


79.6 


07/09/05 


52.5 


92.5 


16/09/05 


19.4 


41.0 


21/09/05 


21.0 


47.8 


Mean 


44.2* 


68.4 


<7 


24.6* 


26.5 




6.6* 


6.8 



*These values are computed without taking into account the 
night of the 05/09/05 (see text for further explanations). 



see that the dispersion of the simulated values (t7 m nh-hi g h 
= 24.6 m and o mn h-iow = 26.5 m) are just slitghly higher 
than the observed a bs = 19.9 m. Also the dispersion of the 
surface layer thicknesses is, therefore, reasonably well recon- 
structed by the model. If we take into account the statistical 
error in the high resolution case (a/v-V), we conclude that 
the mesoscale model provides, for this statistical sample, a 
typical surface-layer thickness just ~ 3 m higher than the ob- 
served one. We have therefore an excellent estimated mean, 
almost within the statistical uncertainty. 



3.2 Optical turbulence vertical distribution: 
seeing in the free atmosphere and in the 
whole atmosphere 

Table [3] reproduces the observed and simulated total see- 
ing (etot) and free atmosphere seeing (efa) for each night. 
The median values with the respective standard deviation 
(a) and statistical error {a/y/N) are reported for the low and 
high resolution cases. One can see that the median total see- 
ing calculated with the low resolution mode (eTOT,mnh-iow 
— 2.05 ± 0.21 arcsec) still remains larger than the observed 
one (tTOT.obs = 1-6 ± 0.2 arcsec) as it has been observed in 
lLascaux et all HH). However the median forecasted see- 
ing in the high resolution mode matches almost perfectly 
the observations (etot \mnh-high = 1-7 ± 0.21 arcsec). The 
same can be stated for the median free atmosphere seeing 
reconstructed by the model (£FA,mnh~hi g h = 0.30 ± 0.17 ar- 
sec) that is very well correlated to the observed one {epA,oba 
= 0.30 ± 0.20 arcsec). Figure [3] shows the correlation plots 
for the total and free atmosphere seeing, for each of the 



Table 3. Total seeing £TOP = £[Sm,h t o P ] an< l seem g in the free 
atmosphere EpA=£[/i ll ,/i lop ] calculated for the 15 nights and av- 
eraged in the temporal range 12:00-16:00 UTC. See the text for 
the definition of h s ; and ht op . In the second column are reported 
the observed values, in the third and fourth columns the sim- 
ulated values obtained with high and low horizontal resolution 
respectively. Units in arcsec. 





Obs. 


MESO-NH 


MESO-NH 






HIGH 


LOW 


Date 


£fa/stot 


sfa/stot 


£fa/ stot 




(/i sI =33m) 


(h s ;=44.2m) 


(h sl =62. 4m) 


04/07/05 


0.3 / 1.6 


0.22 / 2.28 


0.35 / 2.13 


07/07/05 


0.2 / 1.5 


0.28 / 1.91 


0.24 / 1.49 


11/07/05 


1.4 / 1.7 


1.61 / 1.81 


1.57 / 2.08 


18/07/05 


0.3 / 2.0 


0.80 / 1.94 


1.98 / 2.20 


21/07/05 


0.7 / 1.1 


0.86 / 1.27 


0.41 / 1.02 


25/07/05 


0.3 / 1.0 


0.25 / 0.85 


0.23 / 2.03 


01/08/05 


0.8 / 1.6 


0.22 / 2.27 


0.22 / 2.05 


08/08/05 


0.5 / 2.3 


0.35 / 1.70 


0.31 / 1.00 


12/08/05 


0.2 / 1.5 


0.23 / 0.99 


0.28 / 2.49 


29/08/05 


2.5 / 3.6 


2.29 / 2.47 


2.12 / 2.78 


02/09/05 


0.9 / 1.9 


1.16 / 1.54 


0.92 / 1.77 


05/09/05 


0.3 / 1.0 


0.30 / 0.52 


0.27 / 0.73 


07/09/05 


1.4 / 2.8 


1.69 / 3.73 


1.88 / 3.69 


16/09/05 


0.2 / 1.5 


0.21 / 1.57 


0.20 / 2.25 


21/09/05 


0.2 / 1.7 


0.26 / 1.63 


0.22 / 0.75 


Median 


0.3 / 1.6 


0.30 / 1.70 


0.31 / 2.05 


17 


0.7 / 0.7 


0.67 / 0.77 


0.66 / 0.81 


cr/VN 


0.2 / 0.2 


0.17 / 0.21 


0.17 / 0.21 



15 investigated winter nights. This figure confirms the good 
agreement between model and observation from a statistical 
point of view. How do the modifications implemented in the 
Meso-NH model change the the vertical distribution of the 
optical turbulence {C% profiles)? Figure [4] shows the vertical 
profiles of the observed and simulated C% profiles obtained 
with the low and high resolution modes. We can observe 
that, similarly to what has been found in lLascaux et al.l 
(2009), near the ground the high resolution model provides a 
sharper decrease in the optical turbulence than the low res- 
olution mode. However, with this new model configuration, 
the shape of the C% profile is better correlated to measure- 
ments in the first 60 m than with the previous configuration. 
We finally observe that, in both cases, there is still space for 
future improvements in some vertical slabs. We expressly 
avoid more complex calibration spec i fic for a particular site 
as that presented in iMasciadri et al.l 1)20041 '). because we in- 
tend to use the same model to discriminate the quality of 
other sites above the internal plateau. 



4 CONCLUSIONS 

In this paper simulations of the C% profiles related to all the 
15 nights for which measurements done at Dome C in winter 
are available, have been compared to measurements. We 
conclude, that, in the present configuration, the Meso-NH 
model with the 'Astro-Meso-NH' package provides excellent 
estimates of the optical turbulence at Dome C, from a 
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Figure 3. Correlation plot between measured and simulated (a) 
total seeing and (b) free atmosphere seeing. Black dots: low reso- 
lution - Red dots: high resolution. The error bars for simulations 
and observations are reported for one point only (and are equal 
to a). Units are in arcsec. 



qualitative as well as quantitative point of view, if used 
in the high resol u tion m ode. All the conclusions achieved 
m ILascaux et all (|2009h remains unchanged. However we 
clearly put in evidence, in this paper, the positive impact 
of the new numerical configuration of the model Meso-NH 
on the forecast of the optical turbulence, especially for the 
surface layer thickness and the total and free atmosphere 
seeing. In synthesis, the Meso-NH model provides: 

(1) a simulated mean surface layer thickness h a i tmn ,h-high 
— 44.2 ± 6.6 m versus the observed h s i, bs = 35.3 ± 5.1 m. 

(2) a simulated median total seeing ETOT,mnh-hi g h — 1-7 
± 0.21 arcsec versus the observed eroT.oba = 1.6 ± 0.2. 

(3) a simulated median free atmosphere seeing 
£FA,mnh-hi g h = 0.30 ± 0.17 arcsec versus the observed 
EFA,obs = 0.30 ± 0.20. 

We are now ready to apply this model to other antarctic 
regions and to explore the optical turbulence features on oth- 
ers potential interesting sites such as Dome A, South Pole, 
or even Ridge A, in order to perform sites inter-comparison 
and identify the best location for astronomical applications. 
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Figure 4. Median profile measured (black line) with mi- 
crothermal sensors mounted on balloons (from Trinquet et al. 
2008) and simulated with the Meso-NH mesoscale model with 
the low-horizontal resolution (blue line) and the high-horizontal 
resolution (red line), computed over the 15 nights. Up: from the 
ground up to 200 m. Down: from the ground up to 20 km. Sim- 
ulations are considered in the temporal range 12:00-16:00 UTC. 
Units in m -2 / 3 . 
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